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We discuss unquenching of the three-gluon vertex via its Dyson-Schwin¬ 
ger equation. We review the role of Furry’s theorem and present first 
results for the quark triangle diagrams using non-perturbatively calculated 
dressing functions for the quark propagator and the quark-gluon vertex. 
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1. Introduction 

The three-gluon vertex was recently intensively investigated in pure 
Yang-Mills theory m- Its quantitative influence on propagators was stud¬ 
ied, for example, in [3], but also for the quark-gluon vertex it plays an 
important role Consequently, the inclusion of unquenching effects 

constitutes a necessary step for future studies relying on the three-gluon 
vertex. Furthermore, the three-gluon vertex also enters in some beyond 
rainbow-ladder truncations of bound state equations, see, for instance, [TJ. 

In the following we will give a status report on ongoing efforts of in¬ 
cluding quark effects into the three-gluon vertex via its Dyson-Schwinger 
equation (DSE). Essential for these efforts is thereby non-perturbative in¬ 
put for the quark propagator and the quark-gluon vertex, which was studied, 
for example, in B El 

2. Three-gluon vertex Dyson-Schwinger equation 

The full DSE for the three-gluon vertex is depicted in Fig. [l] The 
diagrams with quarks are called the quark triangle and the quark swordfish 
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Fig. 1. The full DSE for the three-gluon vertex. Quarks enter via the quark triangle 
diagrams (7th and 8th diagrams on the right-hand side) and the quark swordfish di¬ 
agram (last diagram). Thick/thin blobs denote dressed/bare vertices. All propaga¬ 
tors are dressed. Wiggly/dashed/continuous lines represent gluons/ghosts/quarks. 




Fig. 2. The truncated DSE for the unquenched three-gluon vertex. 


diagrams. To decouple the three-gluon vertex DSE from the equations of 
higher n-point functions (except for the four-gluon vertex) the following 
truncation is employed: All two-loop diagrams and diagrams with non- 
primitively divergent vertices are neglected. This truncation contains then 
all UV leading diagrams. The truncated equation is shown in Fig. [2j The 
only remaining higher n-point function is the four-gluon vertex, for which 
a model is employed. In Yang-Mills theory this truncation scheme turned 
out to be very successful mm- The qualitative appropriateness of the four- 
gluon vertex model was also confirmed a posteriori , see sans] . However, 
within this scheme the quark swordfish diagram is dropped, which might 
yield important contributions da uni. It remains for future investigations 
to include this diagram by using an appropriate model for the two-quark- 
two-gluon vertex. 

The pure Yang-Mills part of the three-gluon vertex DSE was solved 
self-consistently and with full momentum dependence in [3j [4] within the 
truncation described above. One of the major findings was that the tree- 
level structure is dominant and the other three transverse dressing functions 
are subleading [J]. Therefore, we will consider only the tree-level structure 
as well. In the following we will investigate the contributions from the quark 
triangles using results for the quark propagator and the quark-gluon vertex 
obtained elsewhere. 
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Fig. 3. The quark wave function renormalization Zf(p 2 ) and the quark mass func¬ 
tion M(p 2 ) of the quark propagator. 


2.1. Quark propagator and quark-gluon vertex 

For the quark-gluon vertex a 3PI equation is used instead of its DSE. 
This avoids the appearance of the unknown quark-gluon four-point function. 
The quark-gluon vertex can be decomposed into eight transverse and four 
longitudinal basis tensors. For Landau gauge, we only need the transverse 
part and we decompose the vertex as follows, see [5] for details: 


8 

^im(p,q,p-q) = ^2di{p,q,p -(!) 

2=1 

The tensors i = 1,4, 6 , 7, are chirally symmetric and those with i — 
2, 3, 5, 8 are chirally anti-symmetric. 

The 3PI equation for the vertex is solved together with the quark prop¬ 
agator DSE in the chiral limit. A gluon propagator based on Bang is used 
as an input. Results for the quark propagator, characterized by the quark 
wave function renormalization Zf(p 2 ) and the quark mass function M(p 2 ), 
are shown in Fig. [3] and selected results for the vertex in Fig. [4j Note that 
two of the largest dressing functions are g 2 and ^ 3 , which are only allowed 
due to the breaking of the chiral symmetry. 


2.2. The quark triangles 

There are two quark triangle diagrams with opposite direction of mo¬ 
mentum flow of the quarks, see Fig. [2j This can be interpreted as quarks 
going round in one diagram and anti-quarks in the other. In QED similar 
diagrams exist for the three-photon vertex. However, according to Furry’s 
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Fig. 4. The dressing functions of the quark-gluon vertex at the symmetric point 
corresponding to the eight transverse tensor structures. The left plot shows the 
results for the chirally symmetric and the right plot for the chirally anti-symmetric 
structures. 

theorem these two diagrams cancel each other. The argument can be gener¬ 
alized to any diagram with a closed fermion loop. The general result is that 
there is no photon vertex with an odd number of legs. The case of QED 
must be adapted to the non-Abelian character of QCD, viz., color has to be 
taken into account. The color part can be factorized and reads Tr T a T b T c 
for one diagram, Tr T a T c T h for the other, where T a are the generators of 
SU(3). The color trace yields 



( 2 ) 


In QED the two diagrams cancel, because they are exactly the same except 
for the signs. For QCD, however, the effect of the opposite signs is different, 
since the color parts of the two diagrams have a symmetric and an anti¬ 
symmetric part. The symmetric parts then cancel, as in QED, while the 
anti-symmetric parts add up. We thus need to calculate only one diagram, 
multiply it by a factor of 2 and take only the anti-symmetric color part. 
The fact that d ahc drops is also in agreement with m where it was shown 
that the three-gluon vertex must be proportional to f abc . Hence, we will 
consider only one quark triangle in the following. 

The quark triangle diagram is a static diagram in the sense that it does 
not depend on the three-gluon vertex. Thus it needs to be calculated only 
once and its contribution can be added to the Yang-Mills part during the 
iteration process. This is advantageous, as there are 64 different combi¬ 
nations of dressing functions for the quark triangle diagram from the two 
dressed quark-gluon vertices with 8 dressings each. As an example we show 
the contribution calculated from the dressing gq, corresponding to the trans¬ 
versely projected tree-level tensor, in Fig. [5j This particular contribution 
turns out to be rather small when compared to the full three-gluon vertex 
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Fig. 5. Left : Three-gluon vertex result for the Yang-Mills part from [3]. Right : 
The quark triangle diagram when only the tree-level dressing of the quark-gluon 
vertex, gi, is considered. 


of pure Yang-Mills theory, see also Fig. [5| Nevertheless, it can have quite 
some impact on the three-gluon vertex due to its zero crossing: Even a small 
change can move the zero crossing. Whether this is indeed the case will be 
investigated in future calculations. In particular the contributions from the 
other quark-gluon vertex dressing functions must be included to arrive at a 
final conclusion. 


3. Conclusions and outlook 

The three-gluon vertex plays an essential role for future DSE calcula¬ 
tions in full QCD. Consequently, the inclusion of quark effects is an im¬ 
portant extension of existing pure Yang-Mills calculations. We have set up 
a straightforward scheme that will allow unquenching of the vertex via its 
DSE. Essential input is thereby provided by non-perturbative results from 
the matter part of QCD. However, back-coupling effects of the three-gluon 
vertex onto the matter part are not included in this scheme. 


Acknowledgments 

Funding by the FWF (Austrian science fund) under Contract P 27380- 
N27 and through the Doctoral Program “Hadrons in Vacuum, Nuclei and 
Stars”, Contract W1203-N16, is gratefully acknowledged. 

REFERENCES 

[1] M. Pelaez, M. Tissier, and N. Wschebor, Phys.Rev. D88 (2013) 
125003, arXiv:1310.2594 [hep-th]. 

[2] A. Aguilar, D. Binosi, D. Ibanez, and J. Papavassiliou, Phys.Rev. 

D89 (2014) 085008, arXiv: 1312.1212 [hep-ph] , 








6 


REFERENCES 


[3] A. Blum, M. Q. Huber, M. Mitter, and L. von Smekal, Phys. Rev. D 
89 (2014) 061703(R), arXiv:1401.0713 [hep-ph]. 

[4] G. Eichmann, R. Williams, R. Alkofer, and M. Vujinovic, Phys.Rev. 
D89 (2014) 105014, arXiv: 1402.1365 [hep-ph], 

[5] A. Windisch, Ph.D.Thesis, University of Graz, 2014. 

[6] M. Hopfer, Ph.D.Thesis, University of Graz, 2014. 

[7] C. S. Fischer and R. Williams, Phys.Rev.Lett. 103 (2009) 122001, 
arXiv:0905.2291 [hep-ph]. 

[8] R. Alkofer, C. S. Fischer, F. J. Llanes-Estrada, and K. Schwenzer, 
Annals Phys. 324 (2009) 106-172, arXiv:0804.3042 [hep-ph] , 

[9] E. Rojas, J. de Melo, B. El-Bennich, O. Oliveira, and T. Frederico, 
JHEP 1310 (2013) 193, arXiv: 1306.3022 [hep-ph], 


[10] R. Williams, arXiv: 1404.2545 [hep-ph], 


[11] A. Aguilar, D. Binosi, D. Ibanez, and J. Papavassiliou, . 

Phys. Rev. 


D90 no. 6, (2014) 065027, 

arXiv:1405.3506 [hep-ph] 



[12] M. Mitter, J. M. Pawlowski, and N. Strodthoff, Phys.Rev. D91 
(2015) 054035, arXiv: 1411.7978 [hep-ph], 


[13] M. Pelez, M. Tissier, and N. Wschebor, arXiv: 1504.05157 
[hep-th]. 


[14] D. Binosi, D. Ibanez, and J. Papavassiliou, JHEP 1409 (2014) 059, 
arXiv:1407.3677 [hep-ph], 


[15] A. K. Cyrol, M. Q. Huber, and L. von Smekal, 

Eur.Phys.J. C75 


no. 3, (2015) 102 

, arXiv:1408.5409 [hep-ph] 



[16] M. Hopfer, private communication, 2014. 


[17] M. Q. Huber and L. von Smekal, JHEP 1304 (2013) 149, 
arXiv:1211.6092 [hep-th], 

[18] M. Q. Huber and L. von Smekal, JHEP 1406 (2014) 015, 
arXiv:1404.3642 [hep-ph]. 

[19] N. Smolyakov, Theor.Math.Phys. 50 (1982) 225-228, 










